FERMILAB-PUB-96/212-T 
hep-ph/960823^ 
August 1996 



Renormalization of 
AB = 2 Transitions in the Static Limit 
Beyond Leading Logarithms 



Gerhard Buchalla 

Theoretical Physics Department 
Fermi National Accelerator Laboratory 
P.O. Box 500, Batavia, PL 60510, U.S.A. 



Abstract 

The renormalization group evolution of AB = 2 transition operators is studied at leading 
order in heavy quark effective theory and at next-to-leading order in QCD. We calculate new 
contributions that were not taken into account in previous work and obtain the complete 
result to next-to-leading logarithmic accuracy. The various scheme and scale dependences 
and their cancellation in physical quantities are discussed explicitly. 



The phenomenon of Bd — Bd mixing is one of the most interesting and important effects in 
flavordynamics. Most notably it is sensitive to the mass and the weak coupling V t d of the top 
quark, which play an important role in testing the standard model picture of quark mixing 
in general and of the unitarity triangle and its implications for CP violation in particular. 
Unlike many other loop-induced FCNC processes, Bd — Bd mixing has been observed and 
its measurement provides some constraint on standard model parameters already today. As 
our knowledge of the top quark mass increases and the experimental result for the mixing 
parameter Xd = Am^/T^ gains in accuracy, B — B mixing has the potential of providing 
very valuable information on \ Vtd\- Such information is crucial for a decisive test of our current 
understanding of flavordynamics. 

Unfortunately the analysis and interpretation of B — B mixing is limited on the theoretical side 
through the need to evaluate the hadronic matrix element of a four-quark operator between 
the B meson states. This computation involves QCD long distance dynamics, which pose 
a notoriously difficult problem not yet solved satisfactorily to date. Ultimately the required 
matrix element should be calculable in the framework of lattice gauge theory. However, since 
the b-quark mass is larger than currently attainable lattice cutoffs, it is difficult to treat B 
mesons with fully dynamical quarks in numerical simulations. In this situation it is useful to 
employ heavy quark effective theory (HQET). This approach allows one to extract the depen- 
dence of the matrix element on the heavy mass m b analytically, leaving only light degrees of 
freedom as dynamical variables in a lattice calculation. 

It is the purpose of this Letter to discuss the relation of the full theory operator to its HQET 
counterparts to leading order in the 1/m expansion and to next-to- leading order in renor- 
malization group improved perturbation theory, thereby generalizing, repectively completing 
results that have already been obtained in the literature |2j 

The effective hamiltonian describing Bd — Bd mixing can be written as 

n* B f 2 = ^M^(V t tV td ) 2 C(xM^ t ^ b )Q + h.c. (1) 

with 

Q = {bd) v _ A {bd) v _ A (2) 

and Xt = m 2 /Mw- From this hamiltonian the Bd — Bd mixing parameter Xd = Am^/T^ can 
be derived with the result 

G 2 M 2 

Xd = T B -^—^\V t d\ 2 C(xt([i t ),fi t ,fi b )(Q(Li b )) (3) 

107T Z 171b 

where 

(Qfa)) = (B\Q(ji b )\B) = \f B m 2 B B B {n b ) (4) 

is the required hadronic matrix element mentioned above, fg is the B-meson decay constant 
in the normalization in which /„. = lSlMeU and Bb is the bag parameter, defined through 

(!)■ 

The effective hamiltonian formalism is based on the operator product expansion, which 
achieves a factorization of short- distance and long-distance contributions to the amplitude 
in (H). The dependence on the top-quark mass and the short-distance QCD effects from the 
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high matching scale ji t = 0(m t ) = O(Mw), at which top and W are integrated out, down 
to a scale fib = 0(m&) are contained in the Wilson coefficient function C (xt(fit) , £*t, £*&) ■ This 
coefficient is calculable in RG improved perturbation theory. It is independent of the scale 
fit to the considered order in a s . On the other hand C does depend on the low scale /i& 
and on the renormalization scheme chosen to subtract divergences in its calculation. Both 
dependences are canceled by the hadronic matrix element in the product C ■ (Q). The impor- 
tant feature of the scale and scheme dependence and their cancellation appears only at 0(a s ) 
and hence can be meaningfully addressed only beyond the leading logarithmic approximation. 
For the coefficient function C(xt(fH), fit) £*&) the next-to-leading order analysis has been per- 
formed and discussed in detail in [||. This calculation resums the leading (a" In™ fit/ fib) and 
next-to-leading (a™ +1 In™ fit/ fib) logarithmic QCD corrections to all orders n — 0, 1, 2, ... in 
perturbation theory and eliminates the /^-dependence to order 0(a s ), which considerably re- 
duces the theoretical uncertainty in the short- distance calculation. In fact the /^-dependence 
is entirely negligible after the inclusion of NLO corrections ||. 

Whereas the coefficient C incorporates the short-distance QCD effects from scale fi t down 
to fib, the matrix element (Q) contains the "long-distance" contributions below fit, = 0{rrib)- 
However, since QCD continues to be perturbative from scales fib down to fi ~ lGeV, the 
corresponding contributions to the matrix element may still be extracted in an analytical way 
using the framework of HQET. Within this formalism the full theory matrix element (Q{fib)) 
may be expanded as follows [|TJ 

(Q(fi b ))=C 1 (fib,fx)(Q^)) + C 2 (fib,fi)(Qs^)) + 0(-) (5) 

m 

where we are working to leading order in 1/m. In this order two operators appear in the 
effective theory. They are given by 

Q = 2Ch (+) d)v- A Ch { ~ ) d) v _ A Q s = 2{h^d)s- P {h { - ) d)s-P (6) 

with (hd)v~A = /?-7 M (l — 7s)rf and (hd)s-p = h(l — 7s)rf- Here the heavy quark fields are 
defined in the usual way as 

h&\ x ) = e ±imvx ^-b(x) (7) 

Q, Q s are the HQET operators relevant for the case of a B — » B transition. The field 
creates a heavy quark, while annihilates a heavy antiquark. Since the effective theory 
field (h^) cannot, unlike the full theory field b in Q (0), at the same time annihilate 
(create) the heavy antiquark (heavy quark), explicit factors of two have to appear in (^). Due 
to the mass of the b quark an additional operator Q s with scalar-pseudoscalar structure is 
generated in the effective theory. 

(||) has the form of an operator product expansion. The short-distance contributions of scales 
between fib and /i are factorized into the coefficients Ci^, while (Q(fi)) and (Q s (/i)) contain the 
contributions below /a. In particular the leading dependence of (Q(fi)) on the 6-quark mass is 
fully separated into the coefficient functions and no reference to rrib is present in the effective 
theory matrix elements, which makes them more readily computable by lattice methods than 

<Q(aO>- 

In a first step the coefficients have to be determined at fi = fib by calculating the matching 
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of the full theory operator Q onto its HQET counterparts Q, Q s to 0(a s ) in perturbation 
theory. After performing a proper factorization of long- and short-distance contributions, the 
coefficient functions are independent of the infrared properties of the amplitude @, in partic- 
ular of the treatment of external lines. Therefore, for the determination of the coefficients, the 
arbitrary external states may be chosen according to calculational convenience. The simplest 
technique is to put the heavy quark lines on shell, to set all other external momenta to zero 
and to regularize both ultraviolet and infrared divergences dimensionally. In this case the 
0(a s ) corrections to the matrix elements of Q, Q s vanish. Using this method, performing 
the necessary renormalizations in the MS scheme with anticommuting 75 (NDR scheme) and 
subtracting evanescent contributions in the usual way @], we obtain 



C^HhUb) = 1 + 



4:71 



(7ff-7 (0) )ln— + 



where 



Ci(w, tt.) = ^-B, (9) 

8^ + W -15 B=n iV-l §t = _ 2[N+1) (u) 

and N is the number of colors. For N = 3 this result coincides with the calculation of jlJ.The 
logarithmic term in @ reflects the 0(a s ) scale dependence of the matrix elements of Q and 
Q. Accordingly its coefficient is given by the difference in the one-loop anomalous dimensions 
of these operators, 7^ and ^1 . At the same time this logarithm expresses the complete 
dependence of the full theory matrix element (Q) on the 6-quark mass, to 0(a s ) and to 
leading order in 1/m, which is thus factorized into the short- distance coefficient function C\. 
Besides the logarithm (|8|) contains a scheme dependent constant. We have written this term 
in the form B — B in order to make the cancellation of scheme dependences, to be discussed 
below, more transparent. Here B is a constant characterizing the scheme dependence of the 
full theory matrix element (Q), given here in the NDR scheme (see also [0]). 

In a next step we want to evolve the coefficients from the 6-scale fib = 0(rrib) down 
to a typical hadronic scale // ~ lGeV in next-to-leading order of RG improved perturbation 
theory. In this way we are resumming to all orders the leading (a™ ln n fib/ f^) and next-to- 
leading (a™ +1 \n n fib/fj,) logarithmic contributions, which are of 0{1) and 0(a s ), respectively, 
in the limit fib ^> fi, where the logarithms become large. 

The evolution of the coefficients C T = (C*i, C2) is described by the RG equation 

C(fib,fu)=l T C(fi b ,fu) (12) 



where the matrix of anomalous dimensions of the operators Q, Q s can be written to two-loop 
accuracy as 

7 = S [ i£» ,S> ) + \m) [ ® ® ) (13) 
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The first row of 7 has already been considered to next-to-leading order in and we will use 
the results obtained in this paper here. The vanishing of the (12)-entry in 7 means that Q 
does not mix into the operator Q s . The second row of 7 will be the subject of the present 
article. 

Using the well known general expressions for the solution of (12) at next-to- leading order ||, 
the structure of 7 given in ( |THj ) and the initial values of C\$ from (^), (^j) we obtain after 
inserting the solution C(fib,fi) into (H) 



(QM) = 



a s (ji) 



a s (Hb 
4n 



(7f 1 ) -7 (0) )ln^ + 5 



~(o) 

721 



47T 



7ll 722 
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Ots{Hb) 


di 


a s (nb) 


d 2 \ 




_ a s (n) _ 




_ Ol a (fi) _ 





B-J 



47T 



g s {^b) 

Air 



g s {^b) 



d 2 



B e (QM) 



Here 



where 
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7ii 



Wo 
UN — 2f 



~(0) 
722 
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0i 



J 
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2Po 



N* - —Nf - 2C F f 



(14) 



(15) 



(16) 
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are the beta-function coefficients with iV (/) the number of colors (flavors). 
Eq. (0) represents the relation between the matrix element of the full theory operator Q, 
normalized at the "high" scale fib, and the matrix elements of the leading HQET operators 
Q and Q s , normalized at a low scale fi ~ lGeV. It is valid to leading order in the l/m& 
expansion and to next-to-leading order in RG improved QCD perturbation theory. This 
leading contribution in the heavy quark expansion corresponds to the static limit for the b- 
quark. Going beyond this approximation requires the consideration of several new operators, 
which arise at the next order in 1/m. These contributions have been studied in in the 
leading logarithmic approximation. 

We further remark that (0) as it stands is valid in the continuum theory. In order to use 
lattice results one still has to perform an O(o. s ) matching of Q to its lattice counterpart. This 
step however does not involve any further renormalization group improvement since by means 
of flliD (Q) is already normalized at the appropriate low scale fi. For Q s , whose coefficient is 
0(a s ), a tree level matching is sufficient. 

Note that due to the particular form of the RG evolution and the fact that C2 is only 0(a s ), 
the anomalous dimension coefficients 7^ and 722 in ( |13"D do not appear in (0) and hence are 
irrelevant at NLO. The one-loop anomalous dimension 7^, describing the leading logarithmic 
evolution of the effective theory operator matrix element (Q), has been first calculated in 
and H with the result quoted in fllOD . The computation of the two- loop anomalous dimension 
was performed in pL who obtains in the NDR scheme 



N-l 

' 12N 



57 / 32\ 
127iV 2 + 143iV + 63 - — + f 8iV 2 - 16JV + — J vr 2 



(28iV + 44)/ 



(17) 
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As (|T4J) shows, the complete NLO expression requires in addition to 7^ and 'j^ the one-loop 



anomalous dimensions 722^ and 72° , which describe the evolution of the operator Q s and the 
mixing of Q s into Q, respectively. We have calculated these quantities and find 

^)- iV + 1 *<«>)_ 3iV 2 -4iV-7 

721 " ~^\T~ 722 " N ( 8) 

Z.W AO) 



The determination of 721 and 722 is straightforward (for general reviews see 10|). We only 



mention that in order to reexpress additional operator structures, arising in the calculation, 
in terms of the basis operators Q, Q s , we made use of the identity (valid in four dimensions 
for v 2 = 1) 

^V(l - 75 )® ^7,(1 - 7s) = 47^(1 - 75) ® 7,(1 - 7s) (19) 

and of the projection property jbh^ = of the heavy quark effective fields (|7]). The 

relation (|TJ) is also necessary for the matching calculation leading to @, @. 

Let us discuss the structure of relation (|HD in more detail. In particular we would like to 



consider the scale and scheme dependences present in (p~4]) and to clarify their cancellation in 
physical quantities. The dependence on /1 in the first factor on the r.h.s. of Ql4j) is canceled by 
the //-dependence of (Q(fj,)). The dependence on fi h of this factor is canceled by the explicit 
In fib term proportional to 7^ . Hence the only scale dependence remaining on the r.h.s., to the 
considered order 0(a s ), is the one ~ a s (fib) , y^ ln//&. This is precisely the scale dependence 
of the full theory matrix element on the l.h.s., which is required to cancel the corresponding 
dependence of the Wilson coefficient in (f|). Similarly the term ~ a s (fib)B represents the 
correct scheme dependence of (Q(fJ>b))> while the scheme dependence of a s (fi)J cancels with 
the one of (Q(fi)) and the difference B — J is scheme independent by itself. Likewise the 
remaining contributions involving B s are scale and scheme independent to the considered 
order. 

The relation (|TJp involves two different renormalization scales, the scale fi at which the HQET 
matrix elements are defined and the scale fib where the transition from full QCD to HQET 
is made. As is evident from our discussion, the precise values of these scales are arbitrary as 
long as fi — 0(lGeV) and fib = 0(rrib), where the hierarchy fi <C fib is assumed. However, 
since the logarithms hifib/fi are not really very large in the case at hand, one might consider 
neglecting higher order resummations of logarithms altogether and content oneself with the 
strict a s -corrections alone. This limit can be obtained by expanding (|TJ]) to first order in a s . 
In this case the expression ( |T3D simplifies to the result used in [[!]] when fib is set equal to 
fi. This approximation is fully consistent. Yet it is useful and interesting to have the more 
complete expression flT3|) at hand, which allows one to quantify the effects of the leading and 
next-to-leading logarithms hifib/fi in this relation. 

Numerically, for fi b = 4.8GeV, fi = lGeV, = 0.3GeV and / = 4 flavors, the coefficient 
of a s (fib) I '(^)B s (Q s (fi)) in ( pf) reads 1.13, instead of 1 if the summation of logarithms was 
neglected. The new contribution proportional to B s to the coefficient of (Q(fj)) is about +0.01. 
Although these effects are not very large, they have to be taken into account for consistency 
in a complete analysis within the next-to-leading logarithmic approximation. 



Note added: After completion of this work we became aware of reference [11] in which 



the same topic is addressed. The results of JTTJ] are in agreement with our findings. 
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